home *** CD-ROM | disk | FTP | other *** search
/ c't freeware shareware 1997 / CT_SW_97.ISO / pc / software / wissen / macos / fft121.hqx / FFTs for RISC 1.21 / rffttest.m < prev    next >
Text File  |  1996-06-27  |  330b  |  21 lines
  1. N=64;
  2. j=sqrt(-1);
  3. a=ones(N,1);
  4. row=0;
  5. for i1=1:2:N;
  6.     a(i1) = sqrt(row+i1-1+.77777);    
  7.     a(i1+1) = (row+i1-1+.22222)*(row+i1-1+.22222) / N - N/2;    
  8. end;
  9. a(1)=N+3;
  10. a(3)=1-N;
  11. A2=fft(a);
  12. A=A2(1:N/2);
  13. A(1)=A(1)+j*A2(N/2+1);
  14. a1=reshape(a,N/2,2);
  15. if N<=32
  16.     A
  17.     a1
  18. end
  19. max(abs(A-atrans))./(max(abs(A))/N)
  20. max(abs(a-aitrans))./(max(abs(a))/N)
  21.